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An effective bosonic Hamiltonian of Is excitons with "spin" degrees of freedom in two dimension 
is obtained through a projection procedure, starting from a conventional electron-hole Hamiltonian 
TLeh- We first demonstrate that a straightforward transformation of Ti e h into a Hamiltonian of 
bosonic excitons does not give the two-body interaction between an "up-spin" exciton and a "down- 
spin" exciton, which are created by the left- and right-circularly polarized light beams, respectively. 
We then show that this interaction is generated through a projection procedure onto the subspace 
spanned by Is excitons, as a renormalization effect coming from higher exciton states. The projection 
also renormalizes the interaction between Is excitons with the same spins by a large amount. These 
renormalization effects are crucial for the polarization dependence of the optical responses from 
semiconductors. The present theory gives the microscopic foundation of the phenomenology that 
was successfully applied to the analysis of four-wave mixing experiments in GaAs quantum wells 
strongly coupled to the radiation field in a high-Q micro cavity. 

PACS numbers: 71.35.-y, 71.10.-w, 78.66.-w 



I. INTRODUCTION 



Various theoretical methods have been developed to study optical properties of semiconductors§~@. n T.hese methods 
can be divided into two groups; a "fermionic method" and a "bosonic method" . The fermionic methodtrcl is formulated 
on the Hilbert space of two fermionic species, i.e., photo-excited electrons and holes in semiconductors. In this method, 
one solves the coupled equations of motions, the semiconductor Bloch equations (SBEp, for the particle densities of 
electrons and of holes, and for the expectation value of the polarization of the system. Since this method basically 
relics on the Hartree-Fock (HF) approximation for electrons and holes, it is suitable for higher excitation density, 
where Coulombic screening effects guarantee that exciton correlation effects become less crucial. In order to extend 
this method to lower excitation density, where exciton correlations become important, the three- and four-particle 
correlations should be taken into account rand the truncation scheme should be improved. 

On the other hand, the bosonic methodtJTj is based on excitons, bound states of an electron and a hole. Regarding 
this elementary excitation as a bosonic particle, one constructs an effective Hamiltonian of bosonic excitons, from 
which physical quantities, such as linear and non-linear response functions, can be calculated. Here, the effective 
Hamiltonian should be constructed very carefully, as we will show in this paper. The bosonic method is believed to b e 



valid when the optical excitation is weak and when the photon energy is close to the exciton energy (see section |I A| ) , 
because under these conditions main contributions to the optical properties should come from the excitons created 
in the system. With increasing the photo-excitation intensity, the Coulombic force becomes weaker by increased 
screening and/or increased Fermi energy, and the contribution from free carriers becomes more important. Hence, 
the bosonic method is not valid at high excitation intensity. 

A remarkable feature of the light field as a probing tool of materials is that it has the polarization degrees of 
freedom. Recent experimental studies of semiconductor optics make the best use of this fact to reveal more detailed 
properties of excited states. The polarization degrees of freedom of photons induce "spin" degrees of freedom of 
excitations (see below). Experiments, including the four- wave mixing experiment in the time-domain, have revealed 
the crucial roles of the interaction between an exciton (ox._an eh pair) created by the left-circularly polarized light 
and the one created by the right-circularly polarized lightcHZI. However, most of the existing theories could not treat 
the polarization dependence correctly. For example, since the SBE (in its original form) were discussed within HF 
theoryB, the excitations with different polarization degrees of freedom (e.g., left and right) are completely decoupled. 
Hence it is impossible to explain the polarization dependence of the optical response. It is the case with bosonic 
theories. In most cases, the polarization degrees of freedom were not included ,in the theories for simplicity. Although 
it has been asserted that including such degrees of freedom should be triviaB, the excitons with opposite spins are 
completely decoupled. This is because they essentially end up with HF of Is excitons. This means that to discuss 
polarization dependence in bosonic theories is non trivial. 

In this paper, we show that it is not at all trivial to derive the interaction of excitons created by opposite polarized 
light in a bosonic theory. One crucial result is that to go beyond the HF approximation of Is excitons becomes of 
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paramount importance. One of such investigations is seen in Ref.B, where the modification of exciton binding energy 
is discussed. 

Up to now, a strong objection against the bosonic method was that excitons are not bosonic particles. Moreover, 
until recently there were almost no experimental evidence which verifies an effective bosonic theory both qualitatively 
and quantitatively. However, importaiiLexperimental evidence for the validity of the boson picture has been reported 
recently in a two-dimensional systemoE3. The experiment, discussed in detail below, is a non- linear version of normal 
mode coupling in a high-Q micro cavity. This experiment was stimulated by the fact that optical responses from 
semiconductors that are strongly coupled to the photon field have received much attention in recent years. The strong 
coupling is obtained by confining optically active regions in a high-Q micro cavity, which is made possible by the 
development of nano-structure technology. For linear optical responses,-t|he most noticeable phenomenon is a large 
Rabi-splitting, which has been observed in both inorganicEJ and organic£3 semiconductors (insulators). 

Kuwata-Gonokami et al. have performed fourrwave mixing in a GaAs quantum well (QW) that is strongly coupled 
to the radiation field in a high-Q micro cavityuO. They investigated several polarization configurations under the 
condition that the excitation density is very low, and measured the polariton-polariton scattering signals. The 
experimental results wece, successfully reproduced by a phenomenological model, which is called the weakly interacting 
boson model (WIBMpEHil In the W1BM, excitons are treated as interacting bosons. The good agreement with the 
experimental results demonstrated that the bosonic picture is reliable in the lower-excitation regime. The Hamiltonian 
of the WIBM is 

HwiBM = [ Wca L a k<r + ^e&k Aa + #(4 A CT + fo L a k<r) 

k,cr 
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where photons and excitons are described by boson operators, a^ a and 6k a , respectively, with the spin index a = ±. 
This Hamiltonian has three parameters: W is the interaction strength of the.-£xcitons with opposite spin, R is the 
interaction strength of the excitons with the same spin, and v is the filling factoio. The ratio of these three parameters 
was measured as R : W : gv — —3.0 : 0.2 : 1.0. Once these three parameters are fixed, the experimental results for 
all polarization configurations are fitted very well. This iiidicates that a bosonic method is quite reliable in the lower 
excitation density. Nevertheless, existing bosonic theoriesaO cannot give a microscopic foundation of the WIBM. 
The purpose of our paper is to derive a n effe ctive boson Hamiltonian of f s excitons with "spin" degrees of freedom 



under the three conditions cited in Sec. |H A| . Through this derivation, a microscopic foundation of the WIBM, 
especially the interaction terms of excitons with the scattering strength W and R in the WIBM, is obtained. We 
show that effects of exciton states higher than Is are crucial when deriving the effective Hamiltonian of Is excitonsE_3: 
these effects yield the two-body interaction term corresponding to W in Eq. ([!]), and largely modify the strength of the 
interaction term corresponding to R in Eq. (pi). The higher exciton states are taken into account through a projection 



procedure, and, as shown in Sec. IIB, it is crucial to project the whole exciton space onto the Is excitons subspace 
and to re-construct of the interaction of Is excitons. 

In a more general sense, the bos oniq description of fermionic systems in two-dimension is one of the most attractive 
fields in condensed matter theoryr 6 H 17 l. If the bosonic method is useful in two-dimensional semiconductors, it would 
be one of the examples of the success of a two-dimensional "bosonization" . 

The organization of this paper is as follows: In Sec. II, after explaining the necessity of the projection, exciton bosonic 
operators are introduced and exciton spins are defined. In Sec. Ill, the exciton interaction is discussed without the 
projection, which corresponds to the two-dimensional case of the existing bosonic theorya. It is shown that the exciton 
interaction corresponding to the W term in WIBM is not obtained in this approximation. In Sec. IV, which is the main 
part of this paper, the projection procedure is introduced^, resulting in appearance of the W term and modifications 
of the interaction strength of excitons with the same spins. Discussions and remarks are presented in Sec. V. In 
Sec. VI, the results obtained in this paper are summarized. Throughout this paper, the units % = 1, e = 1, eo = 1 are 
used. 



II. MODEL AND STRATEGY 
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A. Conditions for the validity of the effective theory 



The purpose of the present paper is to derive an effective Hamiltonian of Is excitons in a QW. Here, the "effective 
Hamiltonian" means that it (approximately) describes the optical responses of the QW correctly, although it is a 
function of Is exciton operators only. 

Since any effective theory is valid only in some specific physical situations, we first clarify the physical situations or 
conditions under which we construct the effective theory. In compensation for this limitation, the effective theory is 
quite useful; it gives deep insights by describing the physics simply. In contrast, such insights are hardly obtainable 
from straightforward calculations using the electron-hole Hamiltonian. 

We construct an effective Hamiltonian that describes optical responses of semiconductor QWs under the following 
conditions: (i) The excitation is weak (weak excitation regime) so that the mean distance l ex of photo-created (virtual 
and/or real) excitons is much larger than the Bohr radius ao of the Is exciton; 

lex > ao, (2) 

and (ii) all the photon energies (pump, probe, and signal) tkoi are close to the energy E\ s of the Is exciton; 

\E\ a — hu>i \ <C \Esp — hu>i\ (i = pump, probe, signal) (3) 

where Ei v is the energy of 2p excitons, and (iii) the line width Ti s of Is exciton is smaller than the detuning energies; 

Tig < \Ei s — huji\. (i = pump, probe, signal) (4) 



Physical meaning of these conditions will be explained in section V B . 

Under these conditions, nonlinear optical signals would not be strong in general. One must therefore devise exper- 
imental methods for detecting the signals with a high sensitivity. For this purpose, a genius method was proposed 
by Kuwata-Gonokami et aln, in which an optical cavity with a high Q value is utilized. This point will be discussed 



later in Sec. VD 



B. Necessity of projection and renormalization 



The Hamiltonian of an eh system is defined on the eh Hilbert space H e j, that is an fermionic Hilbert space spanned 
by e and h states. As long as photo excitations are concerned, all excited states are charge neutral. We can thus limit 
ourselves in the charge neutral sector of W e h- The effective Hamiltonian of excitons, which is defined on a bosonic 
Hilbert space spanned by the exciton states, should describe the dynamics of the eh system in this charge neutral 
sector of H e /j. Note that there are two (or more) choices for the bosonic Hilbert space: one is the whole exciton space 
H ex that is spanned by all exciton states, whereas the other is its subspace H ls that is spanned by the Is states only. 
The effective Hamiltonian depends on the choice of the bosonic Hilbert space. 



We will first consider in Sec. Ill the effective Hamiltonian defined on H ex . Its interaction part takes the following 
form: 

W /« = -L^£ J2 Vex(c l ;M;{S})bl + ^ iS / k ,^ 2S2 b k ,^ s M, 3 s 3 + ---, (5) 

{S} M k,k',q 

where 6k v s denotes the excitonic boson operators (defined later) , and is the area of the QW. The first term denotes 
the two-body interactions between excitons of various states, and • • ■ denotes three- and more- body interactions. 
Since exact calculations, which take account of all te rms of Eq. (|5|) to infinite order, are impossible, one has to make 



approximations. Because of conditions (i)-(iii) of Sec. II A, it is tempting to take only the two-body interactions of Is 
excitons, H™*, among many terms of Eq. (||); 

n ^ = ^EE Mq;{" = i*};{s})&U qS X'- q sA'sAs3, (6) 

{S} k,k',q 

where 6k s denotes 6k „s with v = LsEl. Unfortunately, however, we will show later that the replacement -}{f uU ~ mt _ > 
7^1™* is a very poor approximation, which cannot explain the experimental results even qualitatively. This originates 
from the complete neglect of effects of higher exciton states v = 2p, 3d, • • •, which, however, play important roles as 
intermediate states. 
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To resolve this difficulty, we will then consider in Sec. IV the effective Hamiltonian defined on Hi s . It is obtained 



by the projection procedure, by which the dynamics in H ex is projected onto the subspace Hi s that is spanned by 
Is excitons only. In general, a projectio n pro cedure generates dissipative terms in the projected dynamics in the 



subspace. Under conditions (i)-(iii) of Sec. [I A, however, we may neglect the dissipative terms. Namely, the dynamics 



in Hi s can be described by an Hamiltonian dynamics, whose Hamiltonian (effective Hamiltonian) is a function of the 
boson operators 6k s for Is excitons only. Its interaction part, 7y{" h_ ™*, consists of two- and more-body interactions 
among Is excitons. For example, the two-body interaction 7i™* takes the following form; 

^ = ^E E ^(H{^s^-,s>^s.' (7) 

{S} k,k',q 

In this two-body interaction, effects of higher exciton states have been (partly) included as "renormalization effects" , 
which have modified (renormalized) the forms and the strengths of Hf^. Therefore, in contrast to the case of W{™ " int , 
it is reasonable to take 7i™* as an approximation to ~ mt . In fact, we will show that the replacement H{" ~ mt 
— » 7^1™' is a good approximation, which agrees with the WIBM and experimental results. 

In short, one must perform the projection onto the subspace Hi s to get a correct effective interaction of Is excitons. 
The projection procedure modifies both the form and strengths of the effective interaction. This is essential to justify 
the WIBM. In what follows, we will derive the effective interaction 7^™' of Is excitons in a QW from the conventional 
interacting electron-hole Hamiltonian. 

C. Model 

We consider the conduction and the heavy-hole bands in a GaAs QW, which has a direct band gap. We start from 
the following conventional form of the electron-hole Hamiltonian Ti.^: 

2 



H eh = J2 J Ax $( x ) (-^T + ^) 4iW + E 



x / dxdx'^lix^lix^V^i-v'^i^x'^iix). (8) 



Here, V{r) denotes the Coulomb potential, which behaves in a QW of width L as V(r) w e 2 /er for |r| > L, where 
e is the static dielectric constant, and V(r) « constant for |r| < L. The calculation is simplified by taking the limit 
L — ► wherever the singularity at r = is irrelevant. In Eq. (Js[) , "ipe(h)(x) is the field operator of an electron (hole), 
i = {e, /i}, z e (/j) = 1 (—1), x = (r e (^), J*(u\)i J dx = Y]jz J d 2 r^, and similarly for i' and x' . The index denotes 
the z-component of the total angular momentum, which is a good quantum number, when the z-axis is taken in the 
direction normal to the QW layers. The is defined as — 1 times J z of the corresponding valence band electron. In 
a GaAs QW, J% — ±3/2 for a heavy hole, and J| = ±1/2h. A photon with J* h = +1( — 1) creates an electron-hole 
pair with J| = —1/2 (+1/2) and = +3/2 (—3/2) to conserve the total angular momentum. 



D. Strategy 

Since all states that are excited by photons are electrically neutral, the discussion is confined to the charge-neutral 
sector. Then the following exciton operator can be definealj: 



, f , 2 12 1 (■ m e r e +m h r h \ 
V* - E^ J d ^ d rs ^ exp j 



x^(r e - r h ){S\J z e , J z h )M?e, Je)M?h, Jh)- (9) 

Here, the plane wave corresponds to the center-of-mass motion of an electron- hole pair, ^(r) is a wavefunction for 
the e-h relative motion, (S\J*, J£) the Clebsch-Gordan (CG) coefficient, f2 the QW area, and M = m e + rrih- In the 
following, the heavy hole condition < m e <C rrih is assumed, and fJ,(f/) = m e (h\/M. 

From the explicit calculation of the commutation relation for these operators, they can be treated as bosonic 



operators when the particle density is very low. This is satisfied under condition (i) of section [I A . When tp v (r) is the 
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wave function with the quantum number i/(= Is, 2p + , 2p_, • • ■) of a hydrogen atom in two-dimensioncl, the operator 
bqvs is identified with the bosonic operator for an exciton with the relative motion index v. Then, the exciton states are 
labelled by indices q, v and S, where q is the momentum of the center-of-mass motion, v denotes the set of quantum 
numbers for the relative motion v, and S denotes combinations of J* and J^, as shown in Eqs. (|Io| ) and (Al). Since 



S is related to the total angular momentum, it is referred to as "spin" index in the following. For transitions from the 
heavy-hole band to the electron band in a GaAs QW, possible changes of the total angular momentum are AJ 2 = ±1 
and ±2. We here take the final states corresponding to AJ 2 = +1, —1, +2, —2 as S = +, — , a, f3, respectively. They 



are related with \Je,J^} as 



/\+)\ 

l-> 




\a) 


H 


V 1/3) ) 





-1/2, +3/2) \ 
■1/2,-3/2) 
1/2, +3/2) ' 
1/2,-3/2)/ 



(10) 



which should be compared with Eq. (Al) of semiconductors with different band structures. Since the dipole transition 
is associated with AJ 2 = ±1, |+) and |— ) are dipole active, coupling to circularly polarized light with J 2 ^ = ±1, 
whereas |a) and |/3) are dipole inactive. The general form of Eq. (fill) for elliptically polarized light is discussed in 
RefM 

III. INTERACTION OF EXCITONS BEFORE PROJECTION 

In this section, the interaction Hamiltonian corresponding to Eq. (^|) is obtained in order to clarify the difference 
from the two-body interaction obtained through the projection procedure which is discussed in the next section. For 
this purpose, we calculate the scattering amplitude of excitons without any intermediate states. 

A. Is exciton scattering amplitude 

In this subsection, scattering processes which involve only Is excitons are considered without the projection proce- 
dure, resulting to the interaction Hamiltonian of Is excitons through a straightforward transformation. 

Such scattering processes are schematically shown in Fig. 1, where the index "ex.z" should be read as the set of 
indices {k^, z/j, Si}. In this subsection, v L = Is for any i and the index is dropped if any confusion is not expected. 
These processes are composed of two parts: one is a direct process, Fig. 1 (a)~(d) and the other is a fermionic exchange 
process, Fig. 1 (e)~(h). The form of the interaction Hamiltonian of Is excitons is 

H u = ^ E u & {sM^sA'-^sA'^St, (11) 

where the scattering amplitude is written as 

C/(q; {S}) = UZ(q)Ub(Si, S 2 ; S 3 , S 4 ) + U^U^, S 2 ; S 3 , S 4 ). (12) 

The U^ l (q)Ul ) (Si, S 2 ' : S 3 , S4) and U Ex (q)U Ex (Si, S 2 ; S 3 , S4) are the direct and the exchange scattering amplitudes, 
respectively. The expressions of each component in Eq. jl2] ) are 

U£(q)=n J dr e drgdr /l dr^ k+q (r e ,r/ l )^ k /_ q (rg,r' ft ) 

x {V(r e - v' e ) + V(r h - v' h ) - V(r e - r' h ) - V(r' e - r h )} k (r e , r„)0 k , (r' e , v' h ), (13) 
U D (S 1> S 2 ;S 3 ,S 4 ) = ]T (S^J^J^iS^J^J^iSslJ^J^iS^J!',^'), (14) 

^Ex(q) = -° J dr e dr' e dr/ l dr^ k+q (r e ,r h )^ k /_ q (r' e ,r / /t ) 

x {V(r e - r' e ) + V(r h - v' h ) - V(r e - v' h ) - V(r' e - r h )} <£ k (r' e , r A )^(r ei v' h ), (15) 

and 

Ue x (Si, S2; S3, Si) = (Sl\Je 7 Jh){S 2 \Je 7 Jfi ){S?>\ J ? e > Jh) J% , Jh )) (16) 

e ' e ' h ' h 
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where q (r e , r^) is the product of the wave functions of the center of mass- and the relative motion of excitons, 

^(r e ,r h ) = ^=e^' + ^ip(T e -r h ). 
The wave function of the relative motion ip is the Is- wave function of a hydrogen atom in two-dimension 



(17) 



(18) 



where ao is the exciton Bohr radius. In the following, U E ^ Ex -^ and U^, Ex ^. are referred to as "orbital" and "spin" 
parts, respectively. 

First, the orbital parts are calculated. Fourier representations of orbital parts in direct and exchange scattering 
amplitude are 



P1.P2 



l^(Pl)| 2 |<^(P2) 



+^(pi)^*(pi + q)<^(P2)^*(P2 - q) - 2^(p 1 )|^(p 2 )|V*(pi + q) 
- E ^(q + Pi-P2) |^(pi)| 2 |^(p 2 )| 2 

Pl,P2 L 

+<p(pi)<p*(pi + q)<^(p2)^*(p2 - q) - 2(^(pi)|(^(p 2 )|V*(P2 - q) 

Here, the notations with tilde are defined as 



(19) 

(20) 

(21) 
(22) 



Since the transferred momentum in the exciton scattering processes is fairly small, which da of the order of the 
photon momentum, the direct and the double fermionic exchange interactions are negligible^ and the momentum 
dependence of the exchange interaction can be omitted. This allows the approximation U E (q) ^e^q = 0) = 
and U Ex (q) « U Ex (q = 0) = U Ex . The q dependence of the interaction strength is discussed in Ref.CJ. Under these 
condition, U E and U Ex can be directly obtained from (filJ) and (pH) for q — > 0: 



U°n = 0, 



u° Ex = 2J2 v(p pO [I^(p)I 2 ^(p) W) - I^(p)I 2 I^(p') 



(23) 
(24) 



Equation (E3) reflects the charge neutrality of the system. 



B. Scattering amplitude of higher exciton states than Is 

The calculation of the exciton scattering amplitude which includes excitons with v > Is is quite similar to the 
previous calculation involving excitons with only v — Is. For such scattering processes that satisfy the conservation 
1 iws. the amplitude for scattering from initial state (Vi . vo\ to the final state (za>. i/a) is obtained as 

(25) 



previous calculation involving excitons witn omy v = is. tor t-.. , 

laws, the amplitude for scattering from initial state (1/1,1/2) to the final state (1/3,1/4) is obtained as 

^1 (P)*^2 (P)*^3 (p)^4 (P') 

-tf Ul (p)*<fv 2 (p')*<fiu 3 (p)^ 4 (p') 



Ud{v} — Oj 

p,p' 



(26) 
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under the same assumptions as in the previous calculation. Here, yv(p) is the Fourier transform of the corresponding 
wave function in the real space: 

^(p) = -^|d 2 re l P> I ,(r). (27) 



C. Interaction Hamiltonian of excitons before projection 

From the results in previous two subsections, the bosonic Hamiltonian of excitons H = H° + H mt is obtained when 
it is assumed that the exciton density is low, where HP is the free part of excitons. The general form of Ti. lnt through 
a straightforward transformation is 

TJ° Tjs 

n/int _ ST Ex{v}^ Ex{S} A , f h h fOSt> 

K - 2.^ 20, "k+q^sA'-q^S^k'^Ok^Sa- (2») 

kk'qM{S} 

Using this formula, we express H as 

"H = H° + 7~Ci a + T~t'i s + 'Hothers, (29) 

where Uf s + H' ls = H^ 1 , (see Eq. (|6|)), includes the v = Is operators only, and 'Hothers denotes such remaining 
terms as the interaction between Is and 2p excitons, the interaction between 2p and 2p excitons and other similar 
interactions. The Ti.f s consists only of the operators with S = ±, corresponding to dipole active excitons, whereas 
Ti.'i B consists of terms of S = a and f3 operators, including cross terms with S = ± operators. 

From the explicit calculation of the spin part of the exchange scattering amplitude, Ug x r S y, for all combinations of 

{Si, S2, S3, S4} = { + , — , a, f3}, H.f s and H ls are obtained as 

H t = E E & k+qS^-q A' A*' (30) 

S=± kk'q 



kk'q 



E ( o & ki- J'l'-r, J>\<' J'k 



.S=a,/3 



2 U k+qS u k'-qS"k' S u kS 



- & k+q + 6 k'-qS 5 k'S 6 k++ 5 k+q- 6 L-qS 6 k'S 5 k- 



^k +q+ ^k'_ q _V Q 6 k / 3 + ll - C - 



(31) 



where 6 q s = & q isSi an d the effective interaction strength U = Ug x ^ = ^ ls y. 
For d — 2, and in the limit of L — > 0, Eq. ( gjj ) is evaluated as 

U = 2nal - ^) El » lJXk&g, (32) 

where ao is the exciton Bohr radius and is the binding energy of the two-dimensional excitonEl. Note that the 
estimation of Eq. ( p4| ) in the three-dimensional case, using three-dimensional Is wave function, gives the hard-core 
scattering strength obtained in Ref.H. This-|fact shows that the theory presented in this section corresponds to the 
result obtained by the Usui transformationO. 

The point which should be emphasized in the dipole active part is that there is no interaction terms between the 
exciton with S — + and the exciton with S = — (opposite spin exciton interaction), which is represented in the form 
of cx b* + tf_b-b+, and that there are only the interaction terms between the excitons with S = +(—) (equal spin exciton 
interaction). This is due to vanishing of Ug x ^ S y for the corresponding combination of {S}. Here two ways of thinking 
are possible. The one is that the excitons with S = + and the excitons with S = — 1 will not interact with each 
other. The other is that the absence of such an interaction term is due to the fact that so far we calculated to lowest 
order only. A large number of experiments show that excitons with the opposite spins do interact and that such an 
interaction is crucial, for instance, for four-wave mixing in the time-domain. So it is not appropriate to interpret 
Tii s as the effective Hamiltonian for dipole active Is excitons. The straightforward transformation presented in this 
section is not appropriate for deriving the effective Hamiltonian of Is excitons and the projection procedure discussed 
in the next section is indispensable. 
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IV. PROJECTION PROCEDURE 



In this section, the-interaction Hamiltonian of Is excitons, which corresponds to Eq. (Q), is obtained through the 
projection proceduret^l. This procedure yields the correct two-body interactions of Is excitons in the subspace of Is 
excitons. In particular, the interaction term of excitons with the opposite spins is obtained through the projection, 
whereas such an interaction was not obtained in the previous section. In other words, such an interaction is obtained 
when going beyond the HF approximation of Is excitons alone. Moreover, the projection yields a large renormalization 
of the interaction strength of excitons with the same spins. Among works which have the similar motivation of going 
beyqad HF, modification of the exciton binding energy beyond HF is discussed by considering screening effects in 
RefJ. 

In the previous section, the effects of higher exciton states have been completely excluded in deriving the interaction 
Hamiltonian of Is excitons, Eqs. ( |30| ) and ([h]). In this section, the theory is discussed in the subspace spanned by 
Is excitons, where the effects of higher exciton states are renormalized. Since this scenario is quite similar to the 
derivation of an effective Hamiltonian from the Hubbard model with large on-site Coulomb repulsion through a 
projection procedural, the method used here is referred to as projection. In other words, this is nothing but the real 
part of the second order vertex correction in the filed theory, yielding the energy shift . 

The higher exciton states are taken as the intermediate states in the scattering processes of Is excitons. Schemati- 
cally, the scattering processes of excitons shown in Fig. 2 are considered. Since our purpose is to obtain the effective 
interaction of Is excitons, the relative motion indices of four external lines must be v = Is. As for the intermediate 
states, excitons which are connected with Is excitons by dipole transitions are considered. Then, the lowest energy 
excitons for the intermediate states are 2p± excitons. Note that in two-dimensional system, the p states are doubly 
degenerated. 

The Hamiltonian for the Is — v interaction processes is 

Wout = V^Wls-w 

v 

= EEE ^( ( l) 5 L,k+ q ,s 1 6 L,k^q,s 2 Vk',s 4 ^,k.S3 + h.c. (33) 

v {5} kk'q 

Let Horg be the Hamiltonian which includes the kinetic and the interaction terms of Is and v excitons, and with 
eigenstates and eigenenergies defined by 

W org |$) =£?*|$). (34) 

Consider the following Schrodinger equation: 

(Worg+^out)!*) =E\V). (35) 

Symbolically, the solution of this equation is shown as 

l*> = -p%-|*> 

Hi — /tore 



5» 



*'*°«l*> V *„„,!*>, (36) 



$ E - E <S> E ~ ft^B 

where the projection operator V is defined as 

v = \-Y^m- (37) 

The |^) is rewritten as 



where 



a$ = E-E* ' (39) 
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and 



V 



E-Uo 

In lowest order, this wave function can be approximated as 

1 



-Woutl**)- 



E-U 



-Wout|$). 



(40) 



(41) 



From them, the Schrodinger equation for the projected out Hamiltonian, that is, the effective Hamiltonian is obtained 
as 



(E - 23$) a* = X! a *' (®\ n ° 



1 



E-H, 



(42) 



org 



Here, two successive exchange scattering processes are considered, whose individual scattering amplitudes are 
obtained from Eq. (p6[). The resolvent 1/(23 — Ti. org ) is replaced with the difference of the kinetic energy between 
two v excitons and two Is excitons, because the interaction energy of each exciton can be assumed to be small as 
compared to their kinetic energy. The orbital part of the renormalized scattering amplitude of Is excitons, Ug x , 
which is calculated by second order perturbative calculation, is 



TTOl 

u Ex 



1 y 



lff,(K)| : 



K ^ ls 2 {E v + Ky2M) - 2E U 
n 2 (E v + K 2 /2M) - 2E ls 



^y(p-p' + K) 
pp' 



+2#,(p)#.(p - K)<p v (p)<p v (p') - F la (p)P u (p')<pAp - K)^(p' + K) 



As for the spin part, Ug x , it is obtained by the simple product of two successive spin weights Ug x , 



Uex(Si> S 2 ; S3, S4) — U% x {Sx,S%; S, S')U^ X (S, S'; S3, S4). 



(43) 



(44) 



S,S' 



From the explicit representation of the renormalized scattering amplitude U' = U^ X U^ X , the renormalized Hamiltonian 
of Is excitons is written in the form 



His = n ls + n ls + n ls , 



(45) 



where 1i\ s is the Hamiltonian of free Is excitons, and Hf s +H' ls = WQ* , (see Eq. (]?])), includes the v = Is operators 
only. The Hf s consists only of operators with S = ±, whereas H! ls consists of terms of S = a and f3 operators, 
including cross terms with S 1 = ± operators. 

The interaction Hamiltonian which has only dipole active Is excitons Hfg, is 



± U-U' 



20 

w 



y. y. \+as\'-as\> s\s 



S=± kk'q 



~cT y 6 k+ q + & k<- q -v- fo k +! 



(46) 



kk'q 



where U' — is a positive constant which arises from the renormalization of higher exciton states (y = 

2p + , 2p_, •••). Comparing the right-hand side of Eq. ( (46| ) with that of Eq. (|30|), we see that the coefficient of 
the first term is renormalized as U — ► U — U', and that the second term is generated through the projection pro- 
cedure, which leads to the opposite spin exciton interaction. This is due to the fact that though the spin weight 
U% x {+, — ; +, — ) vanishes, the spin weight U% x {+, —,;+,—) 7^ because the intermediate states can make use of the 
state with the 5 = {a, j3}, and U% x {+, — ; a, /3) ^ 0. This shows that the renormalization of higher exciton states 
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results in a renormalized Hamiltonian 7Y ls , which differs, both quantitatively and qualitatively, from the Hamiltonian 
Tifs, where higher exciton states than Is exciton are completely ignored. 

The renormalized interaction Hamiltonian which includes dipole inactive excitons is obtained similarly as 



n' l6 



u-u' 



E E 



2 b i+qS b l-qS b k'sb kS 



U 

H 



S=a,/3kk'q 

+ b l+q + b l -qS b *'sbl.+ +&k + q- & L-qS 6 k'S&k- 

E ( 5 k+qa & k'-q/A'/3frk Q +ll ' C - 
kk'q ^ 



"a E 6 k +q a 6 k'- q /3 6 k'/3^k £ 
kk'q 



(47) 



Since 7i' ls includes dipole inactive Is excitons with S = a, f3, it does not contribute to the optical response in its 
lowest order 



V. DISCUSSIONS AND REMARKS 



A. Microscopic foundation of the WIBM 

Kuwata-Gonokami et aZ.i introduced a phcnomenological Hamiltonian, WIBM, which yields good agreement with 
the experimental four-wave mixing data. The WIBM has two kinds of interaction terms of excitons: One is a repulsive 
term R for the same spin excitons. The other is an attractive interaction W for the opposite spin excitons. As our 
first important result, we note that the phenomenological HamiltonianEi has the same form as our 7^^,, the dipole 
active part of Ti.\ s . This is quite .reasonable because the other part Ti.[ s , which is dipole inactive, should be invisible 
in low-order|-pptical experimentscj. We can, therefore, identify the parameters R and W of the phenomenological 
HamiltonianEi as 

W = (48) 



The value of U', as given by Eq. (|43J), depends on the material parameters such as M and e, and hence is different 
for different materials. It also depends on the QW parameter L. Moreover, when imperfections in the QW are non- 
negligible, the expressions of U' should be modified accordingly. Therefore, even for the same material the rvalues of 
R and W could vary from sample to sample, which seems to be consistent with recent experimental resultsEi Note, 
however, that the existence of both interaction terms of excitons 7if s is independent of such details. 

Since the accurate evaluation of Eq. ( f43| ) is rather tedious, we here estimate the typical value of U' as follows. The 
K-summation in Eq. ( [l3|) is cut off for K > Cl/L (through the cutoff of V) and/or for K > C ao /ao (thtough ip v ), 
where Cl and C ao are cutoff parameters of the order of unity. For the case of the QW sample of Ref. L « ao, 
hence we may cutoff the K-summation for K > C/ao, where C is of the order of unity. For the v summation, we may 
consider v = 2p± states only, because higher exciton states give much smaller overlap integrals. The summations for 
p and p' aro-replaced with integrals and U' is evaluated, as shown in appendix B, as U' w \&.ha^C 2 El® . 

ReferenceQ reported the ratio R : W as 1 : — 15. From Eqs. ( pf8| ) and (|49|), we find that this ratio is reproduced 
by the present theory when the cutoff parameter C ~ 0.3, which is consistent with the requirement that C is of the 
order of unity. Considering that the values of R and W vary sJightly from sample to samplccJ, the agreement seems 
satisfactory. Note that such a small value of R reported in Ref. fa is due to the renormalization ofU — > U—U'. Once the 
agreement of 7if s with, tlicjshenomenological Hamiltonian is established, the agreement with the experiment follows, 
as presented in Refs.&E3E3. That is-iowest-order perturbative calculations for the polariton-polariton scattering 
amplitudes agree with the experimentBEl From the above discussion, we conclude that the present theory yields the 
microscopic background of the two exciton interaction terms in the WIBM. 
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B. Correct effective Hamiltonian of Is excitons 



The correct form of the effective Hamiltonian of Is excitons is the renormalized one instead of the Hamiltonian 
which is obtained without the projection procedure, i.e., 

K f J = nf s + n' ls . (50) 

In facl^the opposite spin exciton interaction in ii.f s , which is absent in H.f s , has been clearly observed experimentally in 
Refsflt]. We have used a low-order perturbation theory to derive Hi s , where successive scatterings in the intermediate 
processes are not considered. However, this does not imply a total neglect of multi scattering processes, because we 
have calculated a Hamiltonian rather than observables. In fact, a systematic summing up of higher order scattering 
processes is incorporated in our theory if one calculates higher-order perturbation terms from H.\ s by writing down 
the Bethe-Salpeter equation. 

Note that Tii s is not positive definite to the fourth order in the exciton operators. The stability of the system should 
be preserved by higher order terms. In general situations, properties of a system described by such a Hamiltonian 
should not be analyzed by a perturbation theory based on the vacuum of the free part. Nevertheless, we can use 
such a perturbation theory in our case, because our exciton theory has the built-in constraint that the ground state 
is the state with no excitons, i.e., the vacuum of 7i\ s . The effective Hamiltonian 7ii s together with this constraint 
constitutes a consistent theory, which justifies the low-order perturbation theory based on the given vacuum, if the 
optical excitation is sufficiently weak. 



C. Comparison with existing bosonic theories 



As mentioned in the introduction, most of the existing bosonic theories treated spinless excitons and argued that 
the effect of spins is obtained by trivial extension of the spinless theoryB. The result of such theories is that excitons 
interact repulsively. The interaction of excitons with opposite spins then is attributed to the biexciton formation. 
However, an explicit expression of an interaction Hamiltonian of Is excitons with opposite spins was not obtainedE3. 
This is quite natural because the conventional method was discussed in the HF approximation of Is excitons only, 
where higher exciton states were completely ignored. Such a theory will not yield the interaction of excitons with 
opposite spins, as shown in this paper. The projection method discussed in this paper proves that one can obtain 
the explicit form of the interaction Hamiltonian of excitons with opposite spins beyond HF approximation of only 
Is excitons. There are several papers where such an approach is motivated!!; however, only the modification of the 
exciton binding energy is discussed. Using our interaction Hamiltonian, one can reproduce these resultscS. 

SeveraLauthors obtain an interaction between S = + and S = — excitons within the HF approximation of only Is 
excitonsEH. In those approaches, exciton operators do not have a spin coordinate, and the exchange interaction is not 
the correct one in the sense that the spin coordinate is not exchanged because of the absence of such a coordinate. 
From this theory, the interaction strength of the excitons with opposite spins is —1 times the one with the identical 
spins. .-The absence of the opposite spin exciton interaction within HF of only Is excitons is confirmed in many 
papersca. 



D. Validity of the theory 



In deriving the effective Hamiltonian of Is excitons, we have assumed three conditions listed in Sec. II A: (i) the 
excitation is weak so that the mean distance of photo-created excitons is much larger than the Bohr radius of the Is 
exciton, and (ii) all the photon energies (pump, probe, and signal) are close to the energy of the Is exciton, and (iii) 
the line width of Is exciton is smaller than the detuning energies. 

The condition (i) allows us to use the boson representation of excitons. This condition may be confirmed experi- 
mentally by the fact that the signal intensity is precisely proportional to the square of the pump intensity, i.e., the 
optical response is well described by ■ This suggests that the two-body scattering processes would be dominant. 
When l ex is increased to l ex ~ ao, the deviation from the boson statistics of the operators & q „s and b„ uS will become 
non-negligible, which invalidates the boson representation^]. 

The condition (ii) means that the Is excitons give dominant contributions. This allows us to project out all states 
higher than Is. 

The final condition (iii) allows us to neglect relaxation process of Is excitons. Namely, the equation of motion of 
the reduced density operator p\ s in the Is exciton subspace generally takes the following form; 
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^ = -[H ls ,p ls ]+fp ls . (51) 
ot l 

Here, 7ii s describes the unitary evolution of p± s , and T is the relaxation operator, which is described by the imaginary 
part of a vertex correction. Since the relaxation processes are less crucial under the condition (iii), one may disregard 
r, i.e., one may consider only the real part of the vertex correction. One might think that under these conditions the 
optical signals would not be strong enough to obtain experimental data. However, Kuwata-Gonokami et al.u proposed 
a genius method to overcome this difficulty: They confined a GaAs QW in a high-Q optical cavity. This results in a 
large splitting of excitonic polariton spectrum. Some nonlinear optical signals are strongest at the polariton energies 
(upper and lower ones) because of the polariton resonance. On the other hand, the dissipation (which creates real 
excitons in the QW) is weak at these energies. [This may be understood by considering possible final states: When the 
initial state is a single photon state (coming from an external light source) that has the energy of the lower-polariton 
peak, the final state cannot be an exciton state (without a photon) because the energy is short to create a real 
exciton.] Therefore, one can obtain strong signals without significant dissipation when the photon energies are close 
to the polariton energies in a high-Q optical cavity. Using this idea, Kuwata-Gonokami et alu measured nonlinear 
optical responses under the conditions (i)-(iii), and showed the validity of the WIBM. Our theory is valid in such a 
case. 

On the other hand, Shirane jgi al. recently demonstrated that the dissipation becomes important when the Q value 
of the optical cavity is loweredEx In this case the relaxation processes of Is exciton becomes important, which means 
that r must be fully considered. These are closely related to the excitation induced relaxation (EID) discussed below. 
The relaxation effect is one of the future problems. 

Another example to which the present theory is applicable may be the optical Stark effect in high-quality samples 
at low temperature. In the past experiments of the optical Stark effect, one had to take the detuning rather large 
(hence, condition (ii) is not satisfied) in order to avoid the absorption tail. The absorption tail would be reduced for 
samples with better quality and at lower temperatures. 



E. Biexciton formation and attractive interaction of excitons 

It has been conjectured! that a "biexciton effect" would be the origin of the U W term", i.e., the opposite spin 
exciton interaction. However, this argument is misleading. The biexciton state is analogous to a hydrogen molecule 
and is formed essentially through the mixing of two Is states having different centers. The mixing of two hydrogen 
atoms yields the bonding and antibonding states, which are represented as (1/ \/2)(c\^c\^ ± ^^c^h^h^, |0). Here, 
c i(2) crea t es an electron in the Is state located at nucleus 1(2), and creates the nucleus. The lower energy state 
is the bonding, that is, molecular state. In the case of excitons with J| = ±1/2 and = ±3/2, the corresponding- 
states are (l/\/2)[&+&! ± b^bj^O), where the k-dependence is omitted in order to focus on the S'-dependence. The 
bonding state (— sign for a positive coupling constant) has a lower energy and is called biexciton. This energy splitting 
between the bonding and antibonding states is induced by an interaction of the form of b\b_b a bf3 + h.c, which is 
included in 7Y' ls (or, before the renormalization, in Ti' ls of Eq. (|3lj)). From a more general point of view, a typical 
bound state in an interacting boson model is the eigenstate 

|bound) = (oW - sign^c^lO), (52) 

of a boson Hamiltonian in the following form: 

U = 0,(0*0 + tfb + c)c + dU) + g{a^tfcd + aV&a). (53) 

The eigenenergy of |bound) is 2u — \g\. Note that the existence of the bound state is independent of the sign of the 
interaction g. 

On the other hand, the W term lowers the energies of both states by the same amount, hence does not play a 
central role in the formation of the biexciton state. The most important effect of the W term is^to lower the energy of 
b\_bl_\0), relative to those of bl_b f + \0) and &1&L|0), and this effect was detected experimentally!! In the framework of 
the present bosonic theory, 7Y' ls lowers the energy of the bonding (biexciton) state relative to that of the antibonding 
state, and thus is crucial for the formation of the biexciton state, whereas the W term lowers the energy of both 
bonding and antibonding states. 
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F. Dipole decoupling and HF approximation 



The original SBE are a theory within HF. There are two reasons that an extension beyond HF approximation 
becomes necessary. One reason is that it is necessary to take exciton-exciton correlation effects into account. When 
the excitation is low, the .optical response from semiconductors are attributed to excitons. Similar features appear 
under high magnetic fields_J, where the most striking feature is that the signal of time-domain four-wave mixing does 
not decay exponentially, whereas the SBE predict a single-exponential decay. These typical two cases are beyond the 
scope of SBE. 

As for the spin degrees of freedom, the excitations created by photons with the opposite circularly polarization are 
completely decoupled within HF. Furthermore, the experimental results of the polarization dependence of four-wave 
mixing signals and quantum beats are not treated within HF. The coupling of the excitation with the opposite spin 
is obtained beyond HF approximation. _____ 

Finally, we discuss the relation .between the fermionic theoriesa__'___ and the present bosonic theory. The HF 
factorization treatment of the SBEcl can not produce the interaction between the excitation created by right-circularly 
polarized light and the excita.tiap-.by left-circularly polarized light. The HF theory, therefore, corresponds to Tlf s , Eq. 
(|30|). It was argued in Rcfs._t___ that the interactions of an exciton with higher states (including free carriers) are 
important, and that the interactions result in the energy shift, the EID, and the "biexcitonic correlations". In the 
bosonic theory in the form of Eq. (p9f), these effects are included in Ti[ s and Ttothers- After the projection is made, 
the relation is roughly as follows. The renormalized Hamiltonian Hf^, Eq. ([hj]), would include the HF term. The EID 
may be described by both T and H' ls . The "biexcitonic correlation" would-be included in H[ s - We belief J. jaat-thc 
present theory thus helps to bridge the gap between the bosonic theories!] and the fermionic theorie»__ of 
e-h systems. However, more detailed comparisons will be a subject of future studies. 

As another future problem, the microscopic-expression of the filling factor v should be discussed on the same footing. 
Although such an expression was derived inl_, it corresponds to the "before projection" in our theory, so that the 
corresponding term "after projection" remains a subject of future research. 



VI. CONCLUSIONS 



In this paper, we have derived the effective Hamiltonian for Is excitons with spin degree of freedom in two-dimension. 
This theory is valid when excitation density is weak and when the photon energy is close to the Is exciton energy, 
because the boson operators of Is excitons are used. Relaxation processes of excitons should be less important like, 
e.g., a QW in micro cavity with high-Q value. What should be most emphasized is that the projection is crucial to 
obtain the effective Hamiltonian of excitons: the correct effective Hamiltonian of Is excitons is not the one which is 
obtained by discarding the all exciton operators of v ^ Is among the full interaction of excitons (see Eqs. (|5|) and (||)), 
because it can not explain the experimental results even qualitatively, as discussed in Sec. III. The higher exciton 
states v = 2p, 3d, ■ ■ ■ play important roles as intermediate states. In order to include such effects, the projection is 
used. Through this procedure, the interactions of excitons with the opposite spins are obtained and the interaction 



strength of excitons with the same spins is drastically modified (renormalized) as shown in Sec. [V. In short, the 
procedure renormalizes both the form and strength of the effective interaction. 

It is also shown that the effective Hamiltonian obtained through the projection provides the microscopic foundation 
of the phenomenological Hamiltonian, Kwibm, proposed inl_. The agreement of the present theory with experiments 
supports the validity of a description of a fermionic system by bosonic fields in two-dimension, if the excitation is 
weak. This is a strong indication that bosonization can be a powerful tool also in higher than one dimension. 

Helpful discussions with Professor Kuwata-Gonokami and Dr. Suzuura are acknowledged. 



APPENDIX A: DIFFERENT BANDSTRUCTURES: THE CASE OF TWO S-TYPE BANDS 

We have obtained the opposite spin exciton interaction within the projection procedure for exciton boson operator 
including electron and hole spin indices. The two Hilbert spaces spanned by the excitons corresponding to left- 
circularly and right-circularly polarized photons are completely orthogonal. 

One might think that this decoupling is due to the spin configuration J| = ±1/2 and — ±3/2 of the electron 
and heavy-hole bands in GaAs. This is, however, not the case: We show that for a different band structure with 
s-type electron and hole bands (Jf = Ji = ±1/2), excitons with £ = ± are also decoupled before the projection. 
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First, we define exciton spins in analogy to Eq. (|10|). Since the J z = J* + — states are twofold degenerate in 
the current case, we consider linear combinations of both states, which are denoted by S = 0, 0'. Now, a spin-matrix 
is defined by 



10') 



1 \ 
10 

l/y/2 1/V2 

V 1/V2 -1/V2 J 



/| + 1/2, +1/2) 
1-1/2,-1/2) 
1 + 1/2,-1/2) 

V 1-1/2, +1/2) 



(Al) 



Note that only the spin variables will be different from the case of = ±3/2, whereas the orbital part remains the 
same. The interaction terms of the dipole active Is excitons with S = ± before the projection procedure are obtained 
as 



rl l/2®l/2 Is 20 



J2 X b l+^s b l-c l s b ^sbks, 



(A2) 



S=± kk'q 



which correspond to Eq. (|30|) . The interaction strength U is shown in Eq. ( |24f ) . Note that the opposite spin exciton 
interaction is absent as in the case of = ±3/2. The interaction term which includes the dipole inactive exciton 
corresponding to Eq. (31) is 



u 



/2®l/2 Is 



kk'q L S=0,0' 



&k'S&k- 
1 



fe k+q- fo L-qS 5 k'S&k- 



6 k+q+ 6 L-q- 6 k'0fek0 - 2 fe k+q + 6 k'-q- & k'0'frkO' + h.C. 
4 6 k+qO & L-qO fo k'0&kO + ^ k +qO' 5 k'-qO' 5 k'0'6kO' + ^Uq & l -q fo k' 0' ^k 0' 
7 & U q 0' b l - q 0' 6 k' O^k + Q^k'-q 0' b * 0' ^k • 



(A3) 



From a calculation analogous to the one in Sec. IV, we obtain the Is exciton interaction terms after projection. The 
projection modifies the dipole active part of the interaction into the form of 

^1/201/2 1" = Of) X X 6 k+qS 5 k'-qS 6 k'S&kS 



S=±kk'q 



~Q X & k+q+V-q 
kk'q 



bl, &k'-^k + , 



(A4) 



which corresponds to Eq. (p6|). Here, U' is the contribution from the higher exciton states, which is the same as 
Eq. d43|). As in the case of Jfc = ±3/2, the opposite spin exciton interaction is present, and the interaction strength 
of the equal spin exciton interaction is reduced. The remaining interaction terms of Is excitons, corresponding to 
Eq. (|47n, are obtained as 



/201/2 Is 



U-U' 



J2 ( fe k+q + & k'-qS 6 k'S&k+ + &l c+q _&L_ qS kk'S&k- 
kk'qS=0,0' ^ 

^)X( 6 k+q+ & k'-q- 6 k'0&kO-6k +q + fcL_ q _&k'0'&kO'+h.C.) 



kk'q 

U-2U' 



An 



X X 6 k+qS fo k'-qS 6 k'S&kS 

S=0,0' kk'q 



TO X ( 5 k+q0 5 k'-q0 6 k'0'&kO' +& k+qO '4'-qO' 5 k'0fckO 



An 

kk'q 

U-U' 

n 



X & k+q0 & k'-q0' & k'0'kk0- 



(A5) 



kk'q 
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APPENDIX B: ESTIMATION OF U' 



The wave functions of 2p states in two-dimension are 

9\/67ra§ 



y(q)2 P ,± 



Qx/y 



Q [l + (3|q|a /2)2]5/2- 



(Bl) 



Assuming that the transferred momentum in each scattering processes is small, the K summation is replaced by the 
products of the number of the states cut-off by C as J2k /(K) ~ (^C 2 / a o) x /(0)< The two arguments in the absolute 
are calculated as 



^n P -p')^( P )^(p)^(p)^(p') = ^(1180 - 819 \og3)^E™, 
p,p' 



^(P ~ P')^L(P)^l S (p')^2p(p)^2p(p / ) = 



127575 7r tij 



3„2 



rp2D 



32768 2£ 



E: 



5 ex 



(B2) 
(B3) 



where a non-dimensional parameter £ is introduced for simplification of the calculation. This parameter is determined 
by the following: The integral shown in the left hand side of the next equation, which is encountered in the above 
calculation, is approximated by the right hand side integral. 



dxJi \ -x) — 
\a ) (1 



( 1+ 2.2)3/2 ( 1 + 9a ,2)5/2 J q 



dxJi ( — x) -p- 
\a J (1 



(l + £x 2 ) 4 



r K ( r_ 
48£ 6 a 3 2 V<y 



(B4) 



where a — ao/2, and pi (x) and Kz{x) are the first order of Bessel function and the second order modified Bessel 
functions, respectively^. The optimum value of the parameter is £ 2 = 5.2. The parameter £ is introduced just for a 
convenience of analytic calculation, and will not change a qualitative nature of the discussion. 
From these, it is rewritten as 



„. i nc 2 n 
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1 



2Eg-2Eg 

~<PU iv)<f>U (P')^2p(p)^2p (P') 



2 ^(P ~ P') ^L(P)^L(P)^2p(p)^2p(p') 

P:P' 

2 



where the first two-times comes from the number of the 2p states in two-dimension. Using E 2 ^ 
E\v = -E 2 ®/9, U' is estimated as U' « 16. 5a^C 2 E 2 ^. 



(B5) 
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FIG. 1. Diagramatically expression of direct (a)-(d) and fermionic exchange interaction (e)-(h) of Is excitons. The holizontal 
lines represent an electron and a hole and the verical line Coulomb interaction of two particles connected by the line. 
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FIG. 2. Successive exchange scattering processes of excitons taken into the projection procedure. The bullet shows the sum 
of the single exchange scattering processes between the Is and v excitons. 
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